Exchange and correlation effects in the transmission phase through a few-electron 

quantum dot 
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The transmission phase through a quantum dot with few electrons shows a complex, non-universal 
behavior. Here we combine configuration-interaction calculations — treating rigorously Coulomb 
interaction — and the Friedel sum rule to provide a rationale for the experimental findings. The 
phase evolution for more than two electrons is found to strongly depend on dot's shape and electron 
density, whereas from one to two the phase never lapses. In the Coulomb (Kondo) regime the phase 
shifts are significant fractions of n (tt/2) for the second and subsequent charge addition if the dot is 
strongly correlated. These results are explained by the proper inclusion in the theory of Coulomb 
interaction, spin, and orbital degrees of freedom. 
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I. INTRODUCTION 

Recent experiments by the Weizmann group have 
stirred much attention to the phase that electrons acquire 
when they traverse a quantum dot (QD) embedded in the 
arm of a Aharonov-Bohm interferometer!^ These fasci- 
nating measurements call for a deeper understanding of 
electron transport through a strongly interacting object. 
If many electrons — say N > 10 — populate a QD, the 
transmission phase O of the tunneling electron displays 
a much studied^ universal behavior: first increases 
by n through the conductance peak and then it lapses in 
the Coulomb valley. Here we focus on the relatively less 
studied few-electron regime, where the phase evolution 
depends on N. In the samples studied in Ref.[T]the phase 
remains constant in the N = 1 valley, independently from 
QD tunings, whereas it lapses in the Coulomb blockaded 
devices of Ref . [2] For N > 2 the phase is not reproducible 
and shows both increments and lapses, likely due to QD 
shape and exchange effects^ 

In spite of several explanations of the few-electron 
scenaricC^HSl a clear picture is still missing. Many 
works i ntrodu ced major simplifications, like spinless 
electrons) 5 * 7 * 8 ! one-dimensional QDs, 9 10 simplified mod- 
els for Coulomb interaction) 5 * 6 * 11 * ^ or other ad hoc 
assumptions!^ Even the interpretation of the simplest 
N = 1 — > ./V — 2 transition is controversial, being vari- 
ably attributed to the occupation of either the same 2 
or a differenlP orbital from that of the first electron, 
to the role of excited doorway channels,™ to electron 
crystallization!^! 

In this paper we compute by fully including exchange 
and correlation effects. The theory is based on the ap- 
plication of the Friedel sum rule (FSR) as generalized in 
Ref. \TE\ to a multiorbital interacting QD — an exact zero- 
temperature result, T = 0. Since the FSR holds for non 
degenerate ground states (GSs) only) 15 * 16 ! it may be ap- 
plied to either singlet Kondo GSs at zero field, B = 0, 
or non-degenerate GSs in the Coulomb blockade regime 
(B > for singlets and B > for doublets and higher- 
spin states). 



It is worth recalling that, in the conductance valleys 
with odd N at T = and B = 0, the QD is always in the 
Kondo regime. In order to reach the Coulomb blockade 
regime by keeping T = 0, one needs to apply the field 
to the Aharonov-Bohm ring to destroy dot-lead Kondo 
correlations, hence removing QD spin degeneracy.^ This 
condition is reached when T -C (J-bB, with fi^ being the 
Bohr magneton and T the QD level width. In the ex- 
periments the temperature is very low (T ~ 30 mK) 
and the field is relatively weak (B ~ 10 mT). There- 
fore, in order to compare measurements with the theo- 
retical predictions reported in this paper, tiny widths T 
are required, which are actually smaller than the values 
presently reported.^ 

As anticipated above, the theory relies on the applica- 
tion of the FSR. The key idea is that the phase variation 
A0 for the addition of one electron to the QD is given by 
integrating the QD spectral density A/"(u>) between two 
consecutive conductance valleys (with AN = 1). This is 
evaluated exactly for an isolated QD via full configura- 
tion interaction (CI) calculations-^! with N < 5. 

The simulations agree with the Coulomb-blockade re- 
sults of Ref. [Tj (i) The phase evolution for 2 < N < 5 
strongly depends on the dot's shape and density, (ii) 
never lapses in the N = 1 valley, independently from dot 
parameters. Additionaly, (iii) in the Coulomb (Kondo) 
regime the increment A0 through the conductance peak 
is significantly smaller than 7r (ir/2) as a consequence of 
strong Coulomb correlation. 

The theory fails to reproduce the Coulomb blockade re- 
sults of Ref. [2j The reason of this discrepancy is presently 
unclear, but it is likely related to the significant depar- 
ture from the condition ruling the range of applicability 
of the theory, i.e., T -C [IbB at T = 0. In order to 
treat the Coulomb blockade regime even at B = 0, one 
could compute the thermal Green function of the fully 
correlated system at T > Tk, with Tr- being the Kondo 
temperature. We leave for future work this alternative 
route, which would immediately provide the transmission 
phase without invoking the FSRp^l 

The structure of this paper is as follows. In Sec. [II] we 
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FIG. 1: (a) Experimental setup, (b) QD single-particle levels 
for a/b < 1. The number sequence points to the consecutive 
filling of six electrons in a non-interacting picture, (c) Trans- 
mission phase and phase shifts A9x<r vs N for the filling 
sequence plotted in (b). 



set the theoretical model and illustrate the usage of the 
FSR to compute the transmission phase variation A® 
between consecutive conductance valleys. After briefly 
recalling the full CI method in Sec. |HI[ we address in 
Sec. |IV| the evolution of the phase in the absence of 
Coulomb interaction (but taking into account the spin 
degree of freedom) . We eventually consider the fully in- 
teracting case in Sec. [VJ 



II. TRANSMISSION PHASE FROM THE 
FRIEDEL SUM RULE 

We model the experimental setup of Ref. 1 as in 
Fig.HTa). Electrons flow along x from source (S) to drain 
(D), tunneling through a QD of elliptical shape. The el- 
lipse is the generic low-energy form for a two-dimensional 
shallow, gate-defined potentialpEl since the lowest-order 
non- vanishing terms of its series expansion are quadratic. 
The scattering matrix S a is diagonal in the spin index 
a in both Kondo and Coulomb blockade regimes. In 
fact, i n the Kondo regime (B = 0), no elastic spin-flip 
occursP^E^l On the other hand, in the Coulomb blockade 
regime only one a channel is active at time for a given 
energy since B removes spin degeneracy. 

Additionaly, we assume^ mirror reflection symmetry 



in the yz plane placed in the QD center, at x = 
[Fig. Ilia)], hence the stationary scattering states ip a {x), 
eigenstates of S a , are either even (e) or odd (o) with 
respect to x — > —x reflection: for |aj| > R, 



Vw(z) = e~ ik ^ + e 



2i9 e<r ik\x\ 



and 



ipoa(x) = sgn(» 



e -»fc|x| _|_ e 2i6 oa e ik\x\ 



Here the even and odd outgoing waves are phase shifted 
by ea and 6* ocr , respectively, k is the wave vector, R is 
the QD nominal longitudinal axis. The left and right 
travelling states il>±j s:rT (x) are superpositions of even and 
odd states: 

Inside the QD, \x\ < R, electrons experience two-body 
Coulomb interactions in addition to the confinement po- 
tential [cf. Eq. @]. 

We first generalize the results of Ref. [20] for spinless, 
non interacting electrons by including spins. The trans- 
mission amplitude t a for travelling states ip±ko-{x) is: 



£(T ^n 



J0 



'8m(e ea -e oa ). (i) 



Each time sin {0 ea — 6 oa ) appearing in Eq. M changes 
sign, due to a variation of either 9 ea or Q oa as a new 
electron tunnels into the QD, then a lapse of n occurs 
for the transmission phase O. Since this happens when 
t a = [cf. Eq. (fTl)], the lapse is located in the conductance 
valleyPUl 

We then include all many-body correlations by con- 
necting the phase shift 6xa per channel (X, a) (X = e, o 
labels the parity) to the exact spectral density Nx<j{u) 
accumulated at the QD via the FSR [cf. Eq. (20) of 
Ref. [H]: 



h dw 



= irjVxa(u)- 



(2) 



Here A/x CT (w) is the density displaced at the QD by the 
electron in the scattering state (X, a) tunneling at the 
energy ftw fixed by the chemical potential /i = fho. We 
mimic the action of the plunger gate in the linear re- 
sponse regime by varying the value of \x with respect to 
the QD energy levels. 

In practice, to use the FSR we integrate it over the en- 
ergy window between two consecutive Coulomb valleys 
with N and N + 1 electrons in the QD, respectively.^ 
Whereas this procedure provides the information on 
the total phase variation only, A9x a = ^ &Nxa = 
■nh f du)J\fxa(w), it allows to compute AMxa- from the 
interacting Hamiltonian of the isolated dot, -Hqd- This 
key result is based on the conservation of the total num- 
ber of scattering plus QD states both in the presence and 
absence of the QD in the arm of the interferometerJiSl 



The CI evaluation of A0x<j relies on the formula 



AAA 



A0 



Xa 



Xa 



ax 



W +1 lOOI 



(3) 



where \^q) is the exact interacting GS of the iso- 
lated QD with N electrons of energy E$ , H QD |*^) = 
Eq \^!q S ), and d a a creates an electron with spin a in 
the orbital of given parity X and further specified by the 
set of quantum numbers ax ■ Equation pi follows from 
Eq. (21) of Ref. H21 which was inferred by connecting 
the phase shift to the delay time spent by the electron 
wave packet in the QD. This delay is obtained by inte- 
grating the wave function square modulus over both time 
and space. By orthogonality of QD orbitals, only terms 
diagonal in ax indices survive in the formula ([3]). 

We eventually link the transmission phase variation 
A0 to AOxcr through Eq. (fTl). In the Coulomb blockade 
regime only one scattering channel (A, a) is active at 
time between two consecutive valleys with respectively 
N and N + 1 electrons. The active channel is univocally 
determined by the total spins and parities of hl 7 ^) and 
|\JJ +1 ), as obtained by CI. On the other hand, in the 
Kondo regime time-reversal invariance (recall that B = 
0) implies that AA/jtt = AA/jfj., evaluated as half the 
Coulomb blockade value gi ven b y Eq. ^. In this way 
we regain at once the resullP^D that A8 = tt/2 for the 
addition of the first electron. In fact, the term on the 
right hand side of Eq. (pi) is trivially one when N = 0. 



III. THE FULL CONFIGURATION 
INTERACTION METHOD 

The interacting Hamiltonian of the isolated QD is 



N 1 



2^ K\r i -r ji 



(4) 



where the single particle (SP) term is 



Pi 
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OyVi) + 7,9*°iHbB. 



H SP (i) 

&rrv a z 

(5) 

Here <Xj = ±1, K is the dielectric constant, m* is the elec- 
tron effective mass, g* is the gyromagnetic factor, and the 
QD confinement frequencies in the x and y directions, u}q x 
and uiQy, have characteristic lengths a — [/i/(m*u;ox)] 
and b = [h/(m*LOQ y )] 1 / 2 , respectively (the ratio a/b is re- 
lated to the ellipse eccentrity). In Eq. Q the weak B 
does not affect orbital degrees of freedom. 

To wholly include in our theory Coulomb correlation, 
we solve numerically the few-body problem of Eq. Q 
by means of the full CI method (also known as exact 
diagonalization, for details see Ref. 18). The CI few- 
body GS [$!$) is essentially a linear combination of the 
Slater determinants N^f/, 



*?>=>>*?> 



with the unknown CjS being the output of the calculation. 
Here the determinants \<&f ) are obtained by filling in all 
possible ways with N electrons the TVgp lowest-energy 
SP orbitals (two- fold spin degenerate at B = 0), eigen- 
states of the SP Hamiltonian (J5J). In the Fock space of 
these Slater determinants i?QD is a large sparse matrix, 
that we exactly diagonalize by means of the parallel code 
DonRodrigo 23 eventually obtaining the coefficients ct of 
Eq. @. 

The diagonalization proceeds in each Hilbert space sec- 
tor labeled by A, the total spin, and the total parity of 
the few-body wave function. After we have obtained the 
GSs |*^) and |*^ +1 ), we evaluate A7V Xct via Eq. (3), 
and eventually infer AO as explained in Sec. [H] 

In the CI calculations reported in Sec. [V] we used 
iVgp = 36 and diagonalized matrices of maximum lin- 
ear size 2.25 xlO 6 . The relative error for the energy was 
less than 1CP 4 for a/b = 1. 



IV. THE SPINFUL NON-INTERACTING CASE 

To illustrate the effect of the inclusion of the spin de- 
gree of freedom in the calculation of the transmission 
phase let us consider for the time being only the 
SP Hamiltonian, -ffsPi and neglect Coulomb interaction. 
The GS is a Slater determinant with the lowest N spin- 
orbitals filled, |*^) = H? =1 cJ^JO) (|0) is the vacuum). 
The Aufbau filling sequence for 1 < N < 6 is depicted 
in Fig. [lib) for a/b < 1. The first two electrons occupy 
the s orbital with opposite spin, then the 3rd and 4th 
electrons fill in the p y orbital, which is shifted in energy 
from the p x orbital by A = fkuo x (l — a 2 /b 2 ). Note that 
the first electron entering a new SP level is always f, 
due to the effect of B. The evaluation of the phase shift 
AOxa at each electron addition is straightforward, since 
in Eq. ([3]) only one addendum gives a non-zero contri- 
bution to A9xa/^ — exactly one — as a new spin-orbital 
(ax, ff) is occupied; the other ones vanish due to the or- 
thogonality of the states. Therefore in Fig. [lib) one has 
AAfxa = 1 for the sequence (A, er) = (e,f), (e,i), (e,t), 
(e,4-)j (<5)t)) (°>4-) of six consecutive electron additions, 
with the p y (p x ) orbital even (odd) under x — > —x. 

The evolution of for the filling sequence of Fig. flTb) 
is shown in Fig.nTc). Both increments and lapses of are 
derived through Eq. (IT]) (lapse locations in the conduc- 
tance valleys with fixed N are arbitrary). A remarkable 
feature of Fig. fTVc) is that increases by n in both tran- 
sitions N = -> N = 1 and N = 1 -*■ N = 2, since 
the first two electrons occupy the same s orbital with 
opposite spin. This is fundamentally different from the 
spinless casepB where a total increase of by 2tt between 
N = and A^ = 2 occurs only if the two electrons occupy 
orbitals of different parities. 

Two lapses of it occur for in the blockaded regions 
with N = 2 and N = 3 [Fig. [ijc)] as the phases et 



(q\ and 9 e ^ increase more than w, respectively. The whole 
pattern of in Fig. fTVc) up to five electrons coincides 




FIG. 2: Contour plot of AO in the (r s , a/b) space, in units of 
7r, for the N — 1 — > N = 2 transition in the Coulomb blockade 
regime. The grey code goes from (white) to 1 (black). The 
value in the Kondo regime is obtained by dividing A0 by two. 



with Fig. 4a of Ref. 1, provided that one interprets the 
N = 3 — > N = 4 smooth transition as a phase lapse (in 
Sec. M we consider an alternative interpretation). This 
agreement is surprising, since in the experiment SP levels 
have a small energy separation (~ 0.5 meV), if compared 
to characteristic Coulomb energies (~ 1 — 3 meV)Pand 
therefore one would expect significant differences from 
the non-interacting model of Fig. [T] On the other hand, 
the O-evolution would be basically the same as in Fig.[T]if 
the interacting QD ground state were well approximated 
by a single Slater determinant, as in Hartee-Fock theory 
where Coulomb interaction is included as a mean field. 



V. THE ROLE OF COULOMB INTERACTION 

When correlation effects beyond the mean-field leveP^ 
are relevant, we expect that AJ\fx a < i, as suggested in 
Ref. [TBI Since the CI ground states \^q) are superposi- 
tions of the Slater determinants |$f ) [cf. Eq. d(5b], after 
expansion on this basis many cross terms give no contri- 
bution to Eq. (pi: the stronger the correlation, the larger 
the number of Slater determinants, the smaller AMxa- 
This seems to be the case in Ref. [T]for A© ~ 37r/4 in the 
N = 1 -)• N = 2 transition of Fig. 4b and AO ~ 3tt/4 
for N — 6 — >• N = 7 in Fig. 5. This could even be 
the case for A6 ~ for N = 3 ->• N = 4 in Fig. 4a of 
Ref. [TJ if one excludes the possibility of a phase lapse. 
Such interpretation is alternative to the one suggested in 
the previous section. 

To assess the impact of correlation in the CI results, 
we parameterize the electron density n of the circular 
dot (wo = uj ax — Lo Qy ) via the dimensionless radius r s 
of the circle whose area is equal to the area per elec- 



tron, r s = l/[a^(wny^ 2 ], where a^ is the effective Bohr 
radius and n is estimated as in Ref. [25l We next fo- 
cus on the evolution of © as a function of both r s and 
ellipse anisotropy ratio a/b. The first electron addition 
in the Coulomb (Kondo) regime always gives A© = tt 
(A9 = 7r/2). Then remains constant in the N = 1 
valley, independently from the values of either r s or a/b 
[cf . Fig. pj] . The second electron addition is analyzed in 
Fig. [2] plotting in the (r s , a/b) space the contour map 
of AG for N = 1 — >• N = 2. Here we vary the anisotropy 
ratio a/b by keeping the ellipse area irab fixed, so the den- 
sity remains constant. The contour lines of Fig.[2]provide 
the value of A® in the Coulomb blockade regime, whereas 
its Kondo counterpart may be simply obtained by divid- 
ing A8 by two (cf. also Fig. [3]). As r s increases, A© 
monotonously decreases, since correlation effects become 
stronger at lower density, as the Coulomb term in Hqd 
[Eq. Q] overcomes the SP term. A similar trend occurs 
by increasing a/b, since a stronger anisotropy effectively 
lowers the dimensionality of the system, again enforcing 
correlation effects P" Note that by overlapping the exper- 
imental value AO ~ 37r/4 with the plot of Fig. [2] we find 
r s m 4 for a/b = 1. This value corresponds to Hujq = 1.2 
meV for GaAs, which is comparable to the experimental 
estimate of 0.5 meV. 1 

From the analysis of CI data of Fig. [2] we find that 
the orbital parities of both the first and second electrons 
are always even, independently from r s and a/b. This 
prediction agrees with the Wigner-Mattis theorem: the 
two-electron GS is always a singlet and the orbital part 
of its wave function is nodelessPZl This result conflicts 
with other explanations pE3J and it is expected to hold 
even in the presence of disorder and/or more complicated 
potentials. 

In Fig. [3] we follow the evolution of O up to five elec- 
trons for a significant range of QD anisotropics. At the 
experimental density (Hluq — 0.5 meV) a slight varia- 
tion of a/b is sufficient to alter the phase behavior for 
N > 2. Indeed, the relative differences between the val- 
ues of a/b for the 3rd, 4th, and 5th panels are as small 
as 5%. Therefore, O is sensitive to fluctuations of the 
experimental QD parameters, as reported in Ref. [TJ 

The occurrence of alternative scenarios in Fig. [3] is an- 
other signature of correlation. In fact, several excited 
states lie very close in energy to the GS, as it is the case 
in the crossover to electron crystallizationPSESl Hence 
a small deformation of the QD shape easily induces a 
crossing between states of different symmetry. We here 
highlight only the most relevant features of a rich zool- 
ogy, focusing on Coulomb blockade results (solid lines in 
Fig. [3]). In a circular QD at such low density (4th panel 
of Fig. [3]) the three-electron GS is a spin quadruplet as 
an effect of correlation PS Because the two-electron GS is 
a singlet, the transition A = 2^A = 3is spin block- 
aded, i.e., AO = without any lapse as ANxa = 0. 
A slight deformation of the QD (3rd and 5th panels) 
changes the N = 3 GS into a doublet, lifting the spin 
blockade (AO ^ between N = 2 and N = 3). The 
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N = 4 GS is a more robust triplet, since the spin po- 
larization is due to Hund's rule — an open shell effect P^ 
However, a stronger deformation of the QD (2nd and 6th 
panels) breaks the orbital degeneracy of the SP levels of 
the 2nd shell inducing a transition to a singlet GS. At 
such anisotropy ratios singlets and doublets typically al- 
ternate for even and odd electron numbers, respectively. 
A further increase of the deformation (fst and 7th panels) 
changes the filling sequence of higher-energy orbitals. 

In Fig. [3] we also plot in the Kondo regime (dashed 
lines) for those cases such that the QD spin is to- 
tally screened by the cloud of opposite-spin tunneling 
electronspH This excludes high-spin GSs other than sin- 
glets an doublets occuring in the 3rd, 4th, and 5th pan- 
els. The hallmark of correlation is that A9 is a frac- 
tion of 7r and 7r/2 in the Coulomb blockade and Kondo 
regimes, respectively [e.g., compare the SP phase evolu- 
tion of Fig. [lTc) with its correlated counterpart in the 
2nd panel of Fig. |3). 

VI. CONCLUSION 

In conclusion, we highlighted the role of exchange and 
correlation in the transmission phase of a few-electron 
quantum dot. Our findings are relevant for transport 
experiments through strongly interacting nano objects, 
including molecules and carbon-based nanostructures. 



FIG. 3: 9 vs IV for different dot anisotropies a/b. Solid 
(dashed) lines refer to the Coulomb blockade (Kondo) regime. 
We used GaAs parameters and huio = 0.5 meV. From top to 
bottom: a/b = (1.5) -1 , (1.15) -1 , (1-05) -1 , 1, 1.05, 1.15, 1.5. 
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